clear
y1=1;
sigma=5;
mu=0.4;
alpha=0.5;
T=3;
g1=0.1;
w2=1;
lambda1=0:0.1:1;
r=length(lambda1);

for i=1:r
    lambda2(i,1)=equas13(mu,sigma,alpha,T,g1,w2,lambda1(i),y1);
    i
end

plot(lambda1,lambda2)

        